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APPLICATION  OF  STATISTICAL  THEORY  TO  BEAM-RIDER 
GUIDANCE  IN  THE  PRESENCE  OF  NOISE. 

II  -  MODIFIED  WIENER  FILTER  THEORY1 
By  Elwood  C.  Stewart 

SUMMARY 


A  study  has  been  made  of  the  application  of  Newton's  modification 
of  the  Wiener  filter  theory  to  the  optimization  of  a  beam-rider  guidance 
system  operating  in  the  presence  of  glint  noise.  Target  and  missile 
motions  are  assumed  to  be  coplanar. 

The  theory  has  been  applied  using  a  typical  variable-incidence 
missile  by  placing  a  realistic  restriction  on  the  mean- square  surface 
deflection  so  that  the  system  operation  is  confined  to  the  linear  range. 
The  transfer  functions  of  the  optimum  guidance  system  are  derived  and 
an  example  missile-control  system  is  synthesized.  It  is  shown  that  the 
minimum  attainable  error  corresponding  to  a  realistic  control  restriction 
is  close  to  that  for  the  Wiener  theory.  Performance  of  the  system  versus 
noise  magnitude  is  given. 

It  is  shown  that  the  most  critical  saturable  quantity  is  control 
deflection  since  constraining  control  deflections  to  realistic  values  can 
prevent  saturation  of  other  important  quantities.  Servo  energy  require¬ 
ments  are  also  greatly  reduced  in  comparison  with  systems  in  which 
saturation  is  allowed  to  occur. 


INTRODUCTION 


Noise  effects  in  missile-guidance  systems  impose  one  of  the  most 
serious  limitations  on  the  effectiveness  of  a  missile.  Noise  signals  in 
general  have  the  effect  of  adding  to  the  system  false  information  which 
cannot  be  distinguished  from  true  information.  Thus  the  missile  responds 
to  the  noise  signals  as  well  as  to  the  tine  signal;  the  miss  distance  is 
thereby  increased.  Since  many  sources  of  noi6e  (such  as  angular  scintil¬ 
lation)  are  inherent  in  the  physical  mode  of  operation  and  cannot  be 
removed,  it  is  important  that  missile  systems  be  designed  to  minimize  the 
miss  even  though  the  noise  is  present.  Statistical  filter  theory  is  useful 
in  performing  this  minimization. 

One  of  the  most  valuable  theories  has  been  developed  by  Wiener 
(refs.  1  and  2).  A  previous  report  has  considered  the  application  of  this 
^•Supersedes  recently  declassified  NACA  Research  Memorandum  A55Ella 
by  Elwood  C.  Stewart.  1955. _ _____ 
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theory  to  the  minimization  of  the  effects  of  radar  glint  noise  in  a  beam- 
rider  guidance  system  (ref.  3)-  It  was  found  there  that  the  optimum  per¬ 
formance  specified  by  the  theory  could  not  be  obtained  because  of  certain 
practical  restrictions,  as  for  example,  limiting  of  control  deflection. 
Further  study  showed,  however,  that  even  in  the  presence  of  limiting,  the 
optimum  performance  could  be  approached  by  the  addition  of  network  filter¬ 
ing  to  the  guidance  system.  On  the  other  hand,  the  optimum  form  for  the 
added  filtering  could  not  be  determined  from  this  study  since  the  system 
desi;$n  wes  based  on  a  theory  which  did  not  take  into  account  the  finite 
range  of  operation  of  the  saturable  quantities.  In  the  present  report  an 
attempt  is  made  to  overcome  this  deficiency  by  determining  the  optimum 
transfer  function  under  this  last  restriction.  The  theory  for  such  an 
approach  was  made  available  in  a  paper  by  G.  C.  Newton  (ref.  4)  on  a 
modification  of  the  Wiener  theory. 

In  the  application  of  this  theory  to  the  missile  guidance  problem  it 
is  necessary  to  make  certain  assumptions.  The  main  assumption  made  in 
reference  3,  that  is,  the  assumption  that  the  target  and  missile  move  in 
the  same  plane,  also  will  be  made  here.  Other  assumptions,  such  as  those 
relating  to  the  class  of  target  maneuvers  and  noise,  are  discussed  in  the 
text. 


SYMBOLS 


Hco 

Hf 

N 

% 

Yo 

a j 
k 

YT 

yn 

Ym 

®id 

®ii 

®T 

€ 


transfer  function  of  optimum  compensating  network 
transfer  function  of  the  fixed  network 

noise  magnitude  or  zero  frequency  spectral  density,  ft2/ radian/sec 
time  constant  of  the  noise  spectrum  shaping  filter,  sec 
optimum  closed-loop  transfer  function 
acceleration  of  target  maneuver,  ft/sec2 

twice  the  average  switching  rate  of  target  acceleration,  l/sec 
target  displacement  from  a  space  reference ,  ft 

apparent  target  displacement  from  true  target  center  due  to  noise, 
ft 

missile  displacement  from  a  space  reference,  ft 

cross-spectral  density  of  the  input  signal  with  the  desired  output 
spectral  density  of  the  input  signal 

spectral  density  of  noise  displacement  yjj,  ft2/radian/sec 
spectral  density  of  target  displacement  y^,  f t2/radian/sec 
error  between  target  and  missile  position,  y^  -  y^j,  ft 
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€<P  component  of  error  e  due  to  target  motion,  ft 

€jj  component  of  error  e  due  to  noise,  ft 

p  Lagrangian  multiplier 

8  control-surface  deflection,  radians 

t  angle  of  yaw,  radians 

p0  optimum  open-loop  transfer  function 

p  open-loop  transfer  function  of  system  approximation  to  p0 

w  angular  frequency,  radians/sec 

(  )  complex  conjugate  of  (  ) 

GENERAL  CONSIDERATIONS 


Glint  noise  is  a  term  that  is  used  to  describe  a  shift  in  the  appar¬ 
ent  target  location  as  determined  by  a  radar.  It  Is  due  basically  to  the 
variable  reflection  characteristics  of  aircraft  targets  and  arises  from 
the  relative  movement  of  the  various  reflecting  surfaces.  Since  the  radar 
utilizes  the  reflected  signal  to  determine  target  location,  variations  in 
the  reflected  signal  are  interpreted  by  the  radar  as  shifts  in  the  target 
location.  Of  the  many  sources  of  noise  which  may  exist  in  a  missile- 
control  system,  glint  noise  is  one  of  the  most  serious  since  it  has  a 
physical  origin  which  cannot  bo  eliminated,  imposing  a  fundamental  limi¬ 
tation  on  missile  effectiveness.  The  situation  in  regard  to  glint  noise 
is  illustrated  in  figure  l(a)  where  the  true  target  position  is  indicated 
as  yT  and  the  glint  noise  is  represented  by  the  displacement  yN. 

The  present  report  is  restricted,  for  the  sake  of  simplicity,  to  a 
two-dimensional  study  in  which  the  target  and  missile  move  in  a  horizon¬ 
tal  plane  as  shown  in  figure  l(a).2  The  attack  situation,  that  is  the 
angle  between  missile  and  target  velocity  vectors,  may  be  arbitrary.  How¬ 
ever,  in  order  to  illustrate  the  theory,  the  tail  chase  illustrated  in 
figure  1(a)  will  be  considered.  It  should  be  noted  that  for  this  case 
the  reference  line  is  fixed  in  space.  The  guidance  system  which  is  con¬ 
sidered  to  be  of  the  beam-rider  type  is  illustrated  in  figures  l(a)  and 
l(b).  Its  function  is  to  make  the  missile  position  coincide  as  closely 
as  possible  with  the  true  target  position  yT.  The  difference  y™  -  yM 
is  indicated  on  this  figure  by  the  error  €,  which  obviously  should  be 
minimized.  The  criterion  of  merit  which  will  be  used  here  is  the  conven¬ 
tional  mean-square  time  average  of  the  error.  This  criterion  is  particu¬ 
larly  appropriate  in  the  case  of  the  beam-rider  system  because  the  target- 
to-missile  range  is  not  normally  transmitted  to  the  missile;  since  the 
missile  never  knows  when  the  target  will  be  reached,  it  is  reasonable  to 
minimize  the  error  for  all  values  of  time. 

2 The  complete  three-dimensional  problem  would  require  a  more  complex 
analysis  than  used  herein.  Possibly  either  the  present  theory  or  Wiener's 
theory  for  multiple  time  series  (ref.  1)  could  be  applied  to  this  case. 
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Emphasis  should  be  made  that  although  the  example  case  is  taken  to 
be  a  tail  chase,  it  contains  the  essential  ingredients  of  the  problem. 
For  other  attack  situations  a  coordinate  system  can  be  defined  such  that 
displacements  and  accelerations  are  measured  from  a  new  reference  line 
which  is  not  fixed  in  space.  This  line  is  one  which  translates  with  con¬ 
stant  velocity  and  without  rotation.  There  are  several  suitable  condi¬ 
tions  under  which  such  a  reference  line  can  be  so  defined.  In  the  first, 
wherein  the  launching  aircraft  flies  a  collision  course  with  the  target, 
the  line  of  sight  between  launcher  and  target  (which  coincides  with  the 
radar  beam,  ideally)  would  be  a  suitable  reference  line.  Another  case 
would  be  one  in  which  the  tracking  radar  is  sufficiently  far  from  the 
target  that  rotations  of  the  line  of  sight  are  small.  In  either  case, 
all  displacement  quantities  as  well  as  accelerations  are  then  measured 
perpendicular  to  the  reference  line  just  as  shown  in  figure  1(a). 


The  design  of  a  system  normally  depends  on  the  inputs  to  be  encoun¬ 
tered,  in  this  case,  the  target  motion  and  the  noise.  Because  of  their 
random  nature,  neither  of  these  quantities  can  be  conveniently  described 
explicitly  as  functions  of  time,  and  statistical  descriptions  are  more 
suitable.  In  the  following  paragraphs  the  target  motion  and  noise  are 
briefly  described. 


The  glint  noise  can  be  defined  by  specifying  both  the  power  spectrum, 
or  power  spectral  density  as  a  function  of  frequency,  and  the  amplitude 
distribution.  Many  measurements  have  been  made  on  the  characteristics  of 
glint  noise.  (For  a  brief  bibliography  see  ref.  3.)  Although  these  meas¬ 
urements  are  somewhat  complicated  and  uncertain,  it  is  generally  found 
that  the  amplitude  distribution  is  approximately  Gaussian  and  that  the 
spectral  density  can  be  adequately  represented  by 


+  1 


(1) 


An  examination  of  glint  noise  spectra  indicates  that  the  break  point 
(1/2«Tn)  is  generally  on  the  order  of  several  cycles  per  second  and  as 
in  reference  3  will  be  taken  to  be  6  cps,  corresponding  to  Tu  =  0.0265 
second.  The  magnitude  of  the  Bpectrum,  N,  of  equation  (l)  depends  on 
factors  such  as  target  size  and  target  aspect  so  that  the  guidance  system 
is  generally  forced  to  operate  over  a  wide  range  of  magnitudes.  This 
range  may  extend  from  7  f t2/radian/sec  for  small  targets  up  to  around 
30  ft2/ radian/sec  for  large  bombers.  For  this  range  of  magnitudes  it 
might  be  thought  that  it  would  be  necessary  to  optimize  the  system  for 
each  noise  magnitude.  As  shown  in  reference  3,  however,  this  wan  not  the 
case  for  the  method  used  therein  since  near  optimum  results  were  obtained 
by  optimizing  only  for  a  mid-range  value  called  the  design  value.  It  is 
not  unreasonable  to  believe  that  the  same  will  hold  for  the  method  of  the 
present  report.  As  in  reference  3,  a  design  value  of  15  ft2/ radian/ sec 
will  be  used  here. 


The  type  of  target  maneuver  upon  which  to  base  the  system  design  can 
never  be  determined  with  certainty,  since  the  target  quite  obviously  may 
maneuver  in  many  different  ways.  A  reasonable  situation  might  be  one  in 
which  the  target  pilot  is  only  aware  of  the  attack  and  therefore  maneuvers 
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in  some  random  manner  to  avoid  toeing  hit.  One  of  the  most  useful  concepts 
in  system  design  and  the  one  which  is  used  in  this  report  is  to  picture  the 
target  evasive  maneuver  as  a  stationary  random  process  in  which  the  target 
turns  at  its  maximum  possible  rate  alternately  in  opposite  directions  with¬ 
out  regard  to  what  the  attacking  missile  is  doing.  The  length  of  time 
between  switches  is  assumed  to  be  randomly  distributed  according  to  a 
Poisson  distribution  [  (l/T)exp(-T/T)] ,  where  T  represents  the  length  of 
the  interval  and  T  the  average  interval  length.  It  can  be  shown  that 
the  spectra.!  density  of  the  target  accleration  is  described  by 

kaip2 

%  =  rt^+k2) 

Here  the  quantity  ay  represents  the  magnitude  of  the  target  accelera¬ 
tion  normal  to  the  beam,  and  k  is  twice  the  average  switching  rate  of 
k  a  2/f.  The  spectral  density  of  the  target  displacement  is  then  given  by 

kay2  .  . 

®T  S  ®*T  =  +  (2) 

For  the  example  tail-chase  situation,  the  target  is  assumed  to  maneuver 
with  ±lg  acceleration  with  an  average  period  of  five  seconds,  which  gives 
ay  =  32.2  ft/sec  and  k  =  0.4  switch/sec. 

There  are  several  important  comments  to  be  made  concerning  this  type 
of  input:  (l)  First  of  all,  a  statistical  description  of  the  target  maneu¬ 
ver  process  is  a  desirable  one,  since  target  motions  cannot  be  described 
as  unique  functions  of  time.  (2)  It  is  clear  that  the  maneuver  assumed  is 
a  severe  one  and  puts  the  system  to  a  good  test;  it  is  often  found  that 
systems  designed  according  to  theories  based  on  either  no  maneuver  or  very 
weak  maneuvers  have  unacceptably  poor  performance  in  the  presence  of  a 
more  severe  maneuver.  (3)  Another  consideration  not  generally  realized  is 
that  the  stationary  process  described  above  is  also  applicable  to  certain 
important  nonstationary  processes.  In  any  real  problem  it  is  apparent 
that  the  inputs  are  distinctly  nonstationary.  For  instance,  they  are  non¬ 
stationary  because  the  target  motion  and  noise  do  not  exist  for  an  infi¬ 
nitely  long  time  into  the  past.  However,  the  nonstationary  character  of 
the  input  is  due  to  the  strict  mathematical  definition.  It  is  clear  that 
in  the  practical  case  it  makes  little  difference  to  the  missile,  so  far  as 
miss  distance  is  concerned,  whether  a  process  persists  over  an  infinite  or 
a  finite  period  so  long  as  the  process  begins  at  a  time  before  the  end  of 
the  attack  by  an  amount  equal  to  or  greater  than  the  system  response  time. 
(Of  course,  the  process  may  terminate  any  time  after  the  attack  is  over 
without  affecting  the  results.)  In  other  words,  an  infinite  period  is, 
for  practical  purposes,  simply  one  which  is  longer  than  the  system  response 
time.  Thus,  when  the  system  response  times  are  short,  results  obtained  by 
means  of  the  stationary  input  apply  directly  to  an  important  class  of  non¬ 
stationary  problems.  The  results  presented  herein  are  in  this  category. 

(4)  It  should  be  pointed  out  that  an  optimum  design  based  on  the  above  tar¬ 
get  maneuver  will  operate  efficiently  against  this  class  of  maneuver  as  a 
whole,  although  not  necessarily  efficiently  against  any  one  particular  tar¬ 
get  maneuver  (e.g. ,  a  single  target  turn).  However,  as  indicated  in  refer¬ 
ence  3,  this  design  will  be  essentially  optimum  for  the  single-turn  maneuver 
as  well. 
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The  primary  objective  of  the  modified  Wiener  theory  is  to  arrive  at 
an  optimum  transfer  function  in  which  the  effects  of  certain  saturating 
elements,  are  considered.  Ideally  it  would  be  desirable  to  determine  this 
optimum  transfer  function  without  any  assumptions  about  whether  or  not  the 
system  would  operate  in  these  saturation  regions.  However,  no  such  theory 
is  available  at  present.  The  approach  considered  herein  is  one  in  which 
the  transfer  functions  of  certain  available  elements  are  chosen  so  that 
saturation  does  not  occur.  Thus  the  system  is  forced  to  remain  linear 
and  linear  methods  are  immediately  applicable.  Although  it  is  not  appar¬ 
ent  at  first  glance,  there  are  certain  indications  that  this  approach  may 
lead  to  desirable  results.  First,  it  appears  intuitively  that  saturation 
is  undesirable  because  of  the  leas  of  both  intelligence  and  ability  to 
control  the  missile  when  in  the  saturated  regions.  Second,  it  is  indi¬ 
cated  in  reference  3  that  the  undesirable  effects  of  saturation  can  be 
reduced  by  the  addition  of  filtering  without  undue  increase  in  the  error; 
this  method  in  effect  tends  to  keep  the  operation  of  the  system  in  the 
linear  range.  These  are  only  indications,  however,  and  the  real  Justi¬ 
fication  in  the  completely  linear  approach  depends  on  the  resultant  error 
performance  which  can  only  be  established  by  investigation. 

The  theory  to  be  used  in  this  approach  follows  Newton's  modification 
of  the  Wiener  filter  theory.  To  describe  the  essentials  of  the  theory  it 
will  be  convenient  to  refer  to  figure  2  where  the  block  diagram  corre¬ 
sponding  to  the  general  filter  problem  is  given.  It  should  be  noted  that 
the  notation  here  corresponds  to  that  in  reference  4.  In  the  general 
situation  there  are  usually  certain  elements  that  are  given  and  may  not 
be  altered  in  the  design  of  the  system.  These  elements  are  denoted  by 
the  transfer  function  in  the  figure.  In  particular  cases  the  fixed 

elements  might  represent  missile  aerodynamics,  servo  motors,  and  so  forth, 
depending  on  the  application.  To  describe  the  limiting  or  saturating 
quantities  it  is  necessary  to  express  them  in  terms  of  the  fixed  elements 
and  either  the  input  or  output  of  the  system.  Which  choice  is  made  de¬ 
pends  on  the  location  of  the  fixed  elements  -  if  the  output  element  is 
fixed,  as  is  usually  the  case,  it  is  desirable  to  express  the  saturating 
quantities  in  terms  of  the  output.  Figure  2  illustrates  a  situation 
wherein  the  output  element  is  fixed.  At  this  point  it  is  necessary  to 
know  the  saturating  quantities.  For  this  purpose  we  introduce  a  ficti¬ 
tious  situation.  Since  in  the  general  case  the  saturating  quantities  are 
different  from  the  input  to  the  fixed  network,  the  input  to  Hj  in 
figure  2  is  imagined  to  be  fed  through  certain  transfer  functions  Hsn 
which  act  as  recording  elements  and  whose  outputs  are  the  saturating 
quantities  0Bn*  It  should  be  pointed  out  that  according  to  this  formu¬ 
lation  of  the  problem  the  outputs  of  the  saturating  elements  do  not  feed 
into  the  system  but  are  used  merely  for  the  purpose  of  evaluation.  From 
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the  above  discussion,  then,  it  is  clear  that  by  specifying  both  the  fixed 
and  saturating  networks,  the  saturating  quantities  can  be  related  to  the 
output.  The  remainder  of  the  system,  indicated  in  figure  2  by  the  compen¬ 
sating  network  Hc0  is  unknown  and  must  be  determined  according  to  a 
specified  criterion  of  merit. 


✓ 


By  way  of  comparison,  the  original  Wiener  theory  considers  only  the 
problem  of  choosing  the  over-all  transfer  function  so  as  to  result  in 
minimum  mean-square  error,  0e,  between  the  actual  output,  90>  and  the 
desired  output,  The  modified  theory,  however,  considers  the  minimi¬ 
zation  of  this  same  mean-square  error  with  a  side  restriction  on  the 
available  range  of  operation  of  the  limiting  quantities;  for  mathematical 
reasons  it  is  convenient  to  consider  this  restriction  to  be  in  the  fora 
of  a  mean-square  limitation.  According  to  the  method  of  Lagrange  (ref.  5) 
this  means  that  the  quantity  to  be  minimized  is  of  the  form 
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where  the  pn's  represent  Lagrangian  multipliers  which  must  be  chosen 
properly  to  restrict  the  operation  of  the  saturable  quantities  08n  to 
very  nearly  the  linear  range. 

In  most  cases  it  is  possible  to  control  saturation  by  an  appropriate 
choice  of  the  mean-square  value  of  the  saturating  quantities  0BD.  In  the 
simple  case  When  a  quantity  has  a  Gaussian  distribution  of  amplitudes  the 
relation  between  the  mean-square  value  of  this  quantity  and  the  probability 
of  its  limiting  is  well  known  and  particularly  simple.  The  distribution  of 
the  saturating  quantity  in  the  actual  caBe  may  not  be  Gaussian  for  either 
of  two  reasons:  First  the  input  may  not  be  Gaussian  or,  second,  there  may 
be  several  saturating  quantities  which  interact.  On  the  other  hand,  even 
in  these  cases  the  distribution  is  frequently  near  enough  to  Gaussian  that 
saturation  can  be  readily  controlled  by  an  appropriate  choice  of  mean- 
square  value. 


Ihe  derivation  of  the  solution  for  the  optimum  compensating  network 
which  minimizes  the  expression  (3)  is  beyond  the  scope  of  this  report, 
but  as  shown  in  reference  4  an  expression  for  this  network  is 


Hf  (ia)4>id(a)e1<rt 
A"  (a) 


da  dt 


where 


>  (4) 


a 

A(u)  -  rHfdtojHfdu)  PnHsndu)HBnd«)l&li(w) 

n*i 


HACA  TN  4278 


7 


the  above  discussion,  then,  it  is  clear  that  by  specifying  both  the  fixed 
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output.  Ihe  remainder  of  the  system,  indicated  in  figure  2  by  the  compen¬ 
sating  network  Hco  is  unknown  and  must  be  determined  according  to  a 
specified  criterion  of  merit. 
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available  range  of  operation  of  the  limiting  quantities;  for  mathematical 
reasons  it  is  convenient  to  consider  this  restriction  to  be  in  the  form 
of  a  mean-square  limitation.  According  to  the  method  of  Lagrange  (ref.  5) 
this  means  that  the  quantity  to  be  minimized  is  of  the  form 
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properly  to  restrict  the  operation  of  the  saturable  quantities  08n  to 
very  nearly  the  linear  range. 


In  most  cases  it  is  possible  to  control  saturation  by  an  appropriate 
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the  saturating  quantity  in  the  actual  case  may  not  be  Gaussian  for  either 
of  two  reasons:  First  the  input  may  not  be  Gaussian  or,  Becond,  there  may 
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in  these  cases  the  distribution  is  frequently  near  enough  to  Gaussian  that 
saturation  can  be  readily  controlled  by  an  appropriate  choice  of  mean- 
square  value. 


The  derivation  of  the  solution  for  the  optimum  compensating  network 
which  minimizes  the  expression  (3)  is  beyond  the  scope  of  this  report, 
but  as  shown  in  reference  4  an  expression  for  this  network  is 
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In  this  equation  the  ®  quantities  represent  certain  combinations  of 
target  motion  and  noise  spectral  densities  (see  Appendix  A),  Hf  refers 
to  the  fixed  network  transfer  function,  Hsn  refers  to  the  transfer  func¬ 
tion  of  the  nth  saturating  component,  and  Pq  represents  the  correspond¬ 
ing  Lagrangian  multiplier.  Of  these  quantities,  all  will  normally  be 
known  with  the  exception  of  the  pQ.  This  means  that  only  the  form  of 
Hco  1b  known  from  this  equation  since  the  numerical  constants  depend  on 
the  values  of  p^.  Hie  pn,  however,  are  implicitly  related  (ref.  4)  to 
the  mean-square  restrictions  on  0sn  by  the  following  equation 


r"s-  r 

“  j' 


Hco  (  iu)HC0  ( ia))Hsn(  iu)H8n(  iu)  *11  (to)  du 


(5) 


Here  It  is  to  be  noted  that  the  integrand  in  equation  (5)  1b  dependent 
on  the  Pi*P2,  •  •  «,Pjj  by  virtue  of  the  dependence  of  Hco  on  these  quan¬ 
tities  (eq.  (4)).  It  can  be  shown  that  in  many  cases  of  interest  the  pn*s 
are  uniquely  determined  by  the  equations  (4)  and  (5)  when  a  suitable  mean- 

square  value,  0sn2,  is  chosen.  This  quantity  is  to  be  chosen  so  as  to 
limit  the  probability  of  saturation  to  a  small  value  according  to  a  cri¬ 
terion  given  later.  In  most  practical  applications  the  complexity  of  the 
operations  involved  in  solving  these  equations  does  not  permit  the  general 
solution  to  be  obtained  explicitly  in  terns  of  the  input  characteristics 
and  fixed  networks.  In  spite  of  this  a  great  deal  can  be  learned  from 
certain  numerical  solutions,  as  will  be  shown. 


Application  of  Modified  Wiener  Theory  to  the 
Beam-Rider  Guidance  System 


Hie  validity  of  the  application  of  the  modified  Wiener  theory  to  the 
beam-rider  guidance  problem  depends  on  Whether  or  not  several  restrictions 
inherent  in  the  theory  are  met.  In  general,  they  are  no  stricter  than  are 
those  of  the  unmodified  Wiener  theory.  First,  the  input  quantities,  target 
motion  and  noise,  must  be  stationary  random  series  (see  ref.  6  for  a 
detailed  definition)  and  defined  by  corresponding  power  spectra  Which  are 
continuous.  It  is  generally  believed  that  displacements  at  the  target  cure 
approximately  stationary  random  series  ( ref .  3)  •  Since  the  beam-rider 
system  operates  from  these  displacements,  the  inputs  to  the  beam-rider 
system  are  also  stationary  random  series.  Second,  the  transfer  function 
of  the  system  must  have  constant  coefficients.  In  general,  however,  the 
kinematic  loop  of  guidance  systems  involves  a  time-variable  range  factor. 

In  particular,  for  the  beam-rider  system  the  time-variable  factor  is  the 
ratio  of  the  launcher-to -missile  and  the  launche  r  -  to  -  targe  t  ranges.  On 
the  other  hand,  because  the  miss  is  determined  primarily  by  What  happens 
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near  the  end  of  flight  during  which  the  variation  in  this  ratio  is  snail, 
it  is  reasonable  to  assume  that  the  requirement  of  constant  coefficients 
is  approximately  met.  Third,  the  transfer  functions  of  the  fixed  and 
saturating  components  must  be  known. 

Sliaplif ications .  -  There  are  many  quantities  which  may  saturate  in  the 
beam-rider  system;  the  most  important  are  (1)  control  motion,  (2)  rate  of 
control  motion,  and  (3)  radar  receiver  voltages .  To  apply  the  theory  to 
all  these  nonlinearities  simultaneously  would  he  a  tremendous  task.  How¬ 
ever,  there  are  indications  that  only  one  limiting  quantity  need  he  con¬ 
sidered  in  the  present  problem. 

One  indication  is  that  of  the  three  types  of  limiting  not  all  are  of 
equal  severity  so  that  one  type  tends  to  predominate.  Intuitively,  it 
might  he  felt  that  control-motion  limiting  is  the  most  serious  type  and 
that  placing  a  restriction  on  the  control  motion  will  satisfactorily  reduce 
the  other  types  of  limiting.  Of  course  the  validity  of  this  presumption 
would  require  verification. 

Another  indication  is  that  of  the  three  types  of  limiting  one  is  more 
fundamental  in  that  it  is  more  difficult  to  remove  than  the  others.  It 
will  he  noted  that  items  (2)  and  (3)  above  are  both  control-system  limita¬ 
tions  while  (l)  is  essentially  an  aerodynamic  or  mechanical  limitation. 
Since  control  systems  are  relatively  versatile  and  can  he  designed  to  cover 
a  wide  range  of  characteristics,  it  appears  that  control-motion  limiting 
is  the  most  fundamental  to  the  problem. 

From  these  indications  it  appears  reasonable  to  consider  only  control- 
motion  limiting.  The  validity  of  the  assumption  that  reducing  this  one 
type  of  limiting  will  reduce  the  other  types  sufficiently  will  then  be 
verified.  In  terms  of  figure  2  this  simplification  means  that  only  one 
saturating-component  transfer  function,  Hs,  is  involved.  If  the  fixed 
network  is  chosen  to  he  the  aerodynamic  transfer  function  so  that  the 
control  motion  is  the  input  to  the  fixed  network,  then  Hs  equals  unity. 
Thus,  a  block  diagram  more  suited  to  the  specific  case  of  the  beam-rider 
system  would  appear  as  in  figure  3.  It  will  be  noted  in  figuie  3  that  the 
input  quantity  called  0i  in  figure  2  has  been  split  into  the  two  parts 
which  exist  in  the  actual  case,  target  motion  yT  and  noise  yN.  Like¬ 
wise  the  desired  output  of  figure  2  is  represented  in  figure  3  by  the 

target  motion  yT.  The  box  HCo>  then,  is  that  portion  of  the  guidance 
system  which  is  to  be  chosen  to  minimize  the  expression  (3)*  It  may  "be 
noted  that  in  the  form  given  the  diagram  does  not  bear  a  direct  resemblance 
to  the  form  of  an  actual  guidance  system  since  the  latter  will  involve 
feedback  loops.  The  form  shown,  however,  typifies  that  of  the  general  fil¬ 
ter  problem  and  is  better  suited  to  calculation.  The  solution  of  this 
problem  can  be  readily  converted  to  that  of  the  actual  beam-rider  system. 
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A  few  additional  remarks  are  appropriate  concerning  the  two  quanti¬ 
ties  and  in  equation  (4).  If  ®n  is  defined  as  the  autocor¬ 
relation  function  of  the  input  in  figure  2,  then 

i  rT 

<Pii(r)  =  lim  /  0i(t)6i(t  +  T)dt 

T+oo  <!i 

In  terms  of  the  input  of  figure  3>  this  "becomes 

i  rT 

91! (t)  =  lim  /  [yT(t)  +  yN(t)][yT(t  +  t)  +  yN(t  +  T)]dt 
T-»oo  -*p 


s  9t(t)  +  cPjj(t)  +  <RrN(T)  + 

It  is  normally  assumed  that  there  is  no  correlation  between  target  motion 
and  noise  so  that  (^(t)  and  9nt(t)  are  zero.  It  follows  that  the  Fourier 

transform  of  9n(T)  is 

$ii(w)  =  $rp(w)  +  $jj(w) 

Since  in  the  present  case  the  target  motion  and  the  desired  output  are 
identical  it  is  easy  to  show  in  a  similar  fashion  that 

$id(w)  = 


General  solution  with  control  motion  restriction.-  By  virtue  of  the 
simplifications  discussed  in  the  previous  paragraphs  the  application  of 
the  theory  to  a  beam-rider-type  guidance  system  becomes  feasible.  For  the 
case  illustrated  by  figure  3,  the  optimum  compensating  network,  Hco,  is 


given  by 


where 


(6) 


In  the  above  equation  the  quantities  4>rp  and  $n  are  known  from  the  tar¬ 
get  motion  and  noise  characteristics  as  previously  discussed.  The  transfer 
function  Hf  of  the  missile  aerodynamics  was  chosen  for  this  study  to 
represent  that  for  a  typical  variable-incidence  missile.  The  significant 
transfer  functions  for  this  missile  are  given  by  the  following  conventional 
equations  (ref.  7) 
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♦  _1_  +  1 
E  *  Td  e(Ta*si  +  2  £aTas  +  l) 

Jfo  =  Td  Vs2  +  2£t)V  +  1 
♦  =  TF  s(Tfcis  +  1) 

from  which  the  fixed  network  becomes 

yM  _  1  t^b2  +  g^b^8  + 1 

f  *  ~6  "  Ts^  sz(TaSs^  +  25aTas  +  l) 


(7) 

(8) 

(9) 


Table  I  summarizes  the  parameters  for  this  particular  missile  at  a  given 
flight  condition.  It  should  be  pointed  out  for  later  comparison  purposes 
that  these  aerodynamics  correspond  to  those  used  in  the  previous  noise 
study,  reference  3  •  The  remaining  parameter  p  is  not  known  so  that 
equation  (6)  can  only  be  used  to  give  the  general  fora  of  HCo«  To 
illustrate  this,  the  fora  of  Hco  is  derived  in  Appendix  A  where  it  is 
shown  that 

2  s2 (Ta282+£  ^aTaS+l) (Ta2s2+2£aTas+l) 

Hco(s)  -  Te  (TpS+1)  (Tr*s2+25rT7s+l)  (Tu2s2+2$uTus+l)  (Tv2s2+2?yTvs+l) 

(10) 

Many  of  the  factors  involved  in  this  equation  are  functions  of  the  unknown 
parameter  p.  The  exceptions  are  (l)  the  aerodynamic  factors,  that  is, 
the  gain  Ts2  and  the  second  quadratic  in  the  numerator,  and  (2)  the  terms 
due  to.  target  motion  and  noise  characteristics,  that  is,  the  first  two 
denominator  factors.  The  remainder  of  the  terms  can  only  be  evaluated 
after  the  proper  value  of  p  has  been  established.  This  value  can  be 
determined  from  the  restriction  placed  on  the  mean-square  control  motion. 
The  analog  of  equation  (5)  in  the  present  case  is 

S2  =  /  Hco(i<‘>)Hco(iu)[®T(<*>)  +  ®n((*))  (ll) 

00 

In  theory,  equations  (10)  and  (ll)  can  be  solved  for  the  compensating 
network,  Hco,  corresponding  to  an  arbitrary  mean-square  control  motion, 

In  thiB  application,  however,  because  of  the  complexity  of  the  fixed 
network  and  the  input  quantities,  the  computation  is  unwieldy.  Alter¬ 
natively  the  following  procedure  will  be  used  here:  First  calculate  the 
function  Hco  from  equation  (10)  for  a  selected  value  of  p.  After  the 
result  is  inserted  in  equation  (ll),  evaluation  of  the  integral  gives  a 
■  t 

value  of  6  .  By  repetition  of  this  process  for  a  few  Judiciously 


12 


NACA  UI  4278 


selected  values  of  p  It  is  possible  to  determine  the  value  of  p 
corresponding  to  the  desired  mean-square  restriction 

As  mentioned  previously,  HCo  does  not  exist,  as  such,  in  the  beam- 
rider  guidance  system.  The  relationship  between  the  filter  problem  being 
solved  and  the  actual  beam-rider  system  is  illustrated  by  a  comparison  of 
figures  3  and  4.  From  figure  3  It  is  apparent  that  if  the  entire  guidance - 
system  transfer  function  is  represented  by  Yo,  then 

Y0  =  HcoHf  (12) 

Thus  the  solution  for  the  compensating  network  HCo  determines  the  desired 
over-all  transfer  function  Yo  for  the  actual  form  of  beam-rider  system 
shown  in  figure  4. 

Unfortunately,  the  theory  does  not  give  the  performance  of  this  opti¬ 
mum  system  directly  in  terms  of  the  resultant  error.  The  minimum  error  can 
be  found,  however,  from  the  following  relationship  (ref.  7) 

e5^  =  eT^  +  ejp 

p  00  poo 

-  /  |1  -  Yo(iu>)|  «t(“>«“  +  /  |Y0(iu)  I  *fldw  (13) 

'•'“00  '-'“CO 

Effects  of  control-deflection  restriction.-  The  procedure  discussed 
in  the  above  paragraphs  has  been  carried  out  for  a  range  of  values  of  the 
Lagranglan  multiplier  p.  The  results  of  this  study  are  presented  in 
figure  5  where  the  errors  due  to  target  motion  and  noise,  the  total  error, 
and  the  control  deflection  are  given  as  functions  of  p.  These  curves 
illustrate  the  nature  of  the  restriction  of  the  control  deflection.  Large 
values  of  p  here  correspond  to  a  large  restriction  and  therefore  small 
control  deflections.  It  is  apparent  from  the  expression  (3)  that  as  p 
decreases  the  restriction  on  control  deflection  is  reduced  until  at  p  =  0 
the  Wiener  case  is  reached.  By  cross-plotting  these  curves  so  as  to  elimi¬ 
nate  the  variable  p  -Which  has  no  significance  in  itself,  the  curve  of 
figure  6  is  obtained.  This  curve  is  fundamental  to  the  problem  in  that  it 
illustrates  the  dependence  of  error  on  available  surface  deflection.  Here 
the  minimum  obtainable  error  according  to  Wiener  theory  is  indicated  by 
the  horizontal  line  off  to  the  right.  The  curve  is  shown  dotted  in  this 
region  in  order  to  indicate  that  these  values  of  error  would  require  unob¬ 
tainable  control  deflections. 

The  importance  of  the  curve  of  figure  6  is  that  it  represents  the  theo¬ 
retical  lower  limit  of  error  corresponding  to  any  restriction  on  root-mean- 
square  control  motion.  Bach  point  on  the  curve  would  be  achieved  by  a  dif¬ 
ferent  guidance -system  transfer  function.  The  interesting  feature  of  this 
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It  is  illuminating  to  trace  to  their  origin  the  various  factors  com¬ 
prising  these  optimum  transfer  functions.  Analysis  shows  that  terms  due 
to  the  basic  Wiener  theory  also  appear  in  the  modified  theory.  In  addi¬ 
tion,  other  terms  due  to  the  aerodynamics  and  the  control-deflection 
restriction  also  appear  in  the  solution.  This  situation  is  illustrated 
by  the  grouping  of  terms  in  equations  (l4),  (15) >  and  (l6).  It  can  be 
said  that  the  terms  •which  arise  from  basic  Wiener  theory  correspond  in 
form  to  the  first  group  of  terms.  However,  certain  constants  are  some¬ 
what  altered,  as  can  be  seen  from  table  II  by  a  comparison  of  the  param¬ 
eters  given  for  the  two  theories  for  identical  target  motion  and  noise 
characteristics.  The  second  group  of  terms  in  the  equations  indicates 
that  the  optimum  form  of  the  filtering  term  due  to  control-motion  restric¬ 
tion  is  quadratic .2  The  last  group  of  terms  in  the  transfer  functions  can 
be  traced  to  the  aerodynamics.  In  fact,  the  numerators  are  identical  with 
terms  in  the  aerodynamic  transfer  function. 

The  transfer  functions  are  illustrated  in  figures  7(a),  7(b),  and 
7(c);  these  curves  represent  the  optimum  transfer  functions  in  the  pres¬ 
ence  of  the  control-motion  restriction.  For  comparison  purposes  the  cor¬ 
responding  transfer  functions  obtained  from  the  Wiener  theory  are  also 
plotted  on  these  figures.  It  can  be  seen  that,  in  general,  the  transfer 
functions  for  these  two  cases  are  similar  at  the  very  low  frequencies  and 
separate  to  a  greater  extent  as  the  frequency  is  increased.  The  difference 
is  primarily  due  to  filtering  terms  which  are  required  to  satisfy  the 
restriction  on  control  motion. 

The  added  complexity  has  certain  effects  on  the  control  motion  and 
error.  These  effects  can  best  be  illustrated  by  means  of  the  correspond¬ 
ing  power  spectra,  that  is,  by  the  integrands  in  equations  (11)  and  (13). 
These  spectra  have  been  plotted  in  figures  8  and  9,  and  for  comparison  the 
corresponding  curves  for  the  Wiener  theory  are  shown.  It  can  be  seen  from 
figure  8  that  the  spectrum  of  the  control  deflection  for  the  more  complex 
system  attenuates  far  more  rapidly  than  does  that  for  the  Wiener  theory. 

It  is  as  a  result  of  this  fact  that  the  modified  system  does  not  limit. 

As  for  the  error,  it  is  seen  from  figure  9  that  over  the  important  fre¬ 
quency  range  the  error  due  to  noise  is  not  altered  too  seriously  so  that 
most  of  the  increase  in  the  total  error  is  due  to  an  increase  in  the  tar¬ 
get  motion  component. 

Effect  of  noise  magnitude  on  performance.-  The  preceding  results  were 
determined  on  the  basis  of  a  design  noise  magnitude  which  lay  between  the 
expected  extremes.  The  noise  magnitude  may  vary  because  of  targets  of 
different  size,  aspect,  and  turning  rate.  For  this  reason  the  error  per¬ 
formance  for  the  optimum  transfer  function  has  been  determined  and  is 

3It  is  interesting  to  note  that  from  calculations  not  presented  here 
it  has  been  observed  that  the  damping  ratios  of  these  quadratics  always 
lie  between  0.7  and  0.8  over  the  entire  range  of  control  motion  presented 
in  figure  6. _ 
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shown  in  figure  10  (curve  A)  as  a  function  of  the  zero  frequency  noise 
magnitude,  N.  The  error  performance  shown  is  optimum  only  at  the  design 
value  of  noise.  Shown  for  comparison  is  curve  B,  which  represents  the 
error  performance  derived  from  the  Wiener  theory,  corresponding  to  the 
same  design  condition  (see  ref.  3).  The  difference  between  these  two 
curves,  then,  is  the  increase  in  the  error  as  a  result  of  the  rms  restric¬ 
tion  on  the  control  deflection. 

As  a  further  comparison  two  additional  curves  obtained  from  refer¬ 
ence  3  are  shown.  First,  curve  C  illustrates  the  performance  of  a  system 
(with  certain  saturable  elements)  which  was  optimized  for  the  fastest  pos¬ 
sible  response  to  a  step  input  in  the  absence  of  noise.  The  performance 
can  be  seen  to  be  significantly  poorer  than  the  optimum  performance  given 
as  curve  A.  Second,  the  best  result  that  could  be  obtained  in  reference  3 
is  shown  as  curve  D  in  the  figure.  The  system  corresponding  to  this  curve 
was  not  optimum  because  the  linear  Wiener  theory  used  in  this  approach  was 
not  capable  of  evaluating  all  of  the  various  filtering  terms  discussed  in 
preceding  paragraphs.  For  example,  terms  due  to  Wiener  theory  were  not 
altered,  and  terms  due  to  aerodynamics  were  not  used.  Terms  due  to  con¬ 
trol  restriction  were  included  but  the  performance  was  limited  by  the 
arbitrarily  assumed  form.  Thus  the  system  was  not  optimum  and  limiting 
occurred;  the  performance  deteriorated  accordingly. 

System  synthesis.-  To  achieve  the  results  which  have  been  presented 
it  is  necessary  to  design  the  guidance  system  shown  in  figure  4  to  have 
the  optimum  transfer  function  Y0.  This  can  be  done  in  many  different 
ways  since  the  actual  beam-rider  system  consists  of  two  distinct  parts: 
the  tracking  radar  and  the  missile-control  system.  For  example,  the  opti¬ 
mum  transfer  function  Y0  might  be  split  up  and  apportioned  between  these 
two  parts  of  the  guidance  system.  Or,  on  the  other  hand,  the  missile- 
control  system  alone  might  be  designed  to  approximate  the  optimum  Y0  in 
which  case  the  tracking  radar  should  be  designed  to  have  a  relatively  fast 
response.  In  a  similar  manner  the  tracking  radar  could  be  optimized,  in 
which  case  the  fast  response  should  be  designed  into  the  missile-control 
system. 

To  investigate  all  these  possibilities  is  beyond  the  intended  scope 
of  this  report.  As  an  example,  however,  of  one  possible  design  the  choice 
of  missile-control-system  filtering  was  investigated.  For  this  case  it  is 
desired  to  design  the  missile-control  system  to  match  the  optimum  transfer 
function  Yo«  ThiB  procedure  is  not  unique;  the  analysis  used  herein  is 
based  primarily  on  cut-and-try  procedures.  It  has  been  found  that  the 
desired  system  could  be  synthesized  in  a  conventional  fashion  illustrated 
in  figure  11.  To  show  how  this  system  is  capable  of  supplying  the  optimum 
filtering  the  open-loop  transfer  function  has  been  derived  in  Appendix  B 
where  it  is  shown  that  with  certain  assumptions 
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?il  o  k2kBv 
5  I'd  + 


(T10e  +  l)(T4s  +  l) (Tfr2s2  4  2gbTbS  +  1) 
82(ThS  +  1)(T58  +  l)(as®  +  bs2  +  cs  +  1) 


(17) 


Comparison  vith  the  optimum  uG  In  equation  (16)  shows  that  certain  dif¬ 
ferences  In  form  exist.  However,  equation  (l6)  can  he  expressed  in  another 
form  by  factoring  approximately  certain  quadratic  terms  which  have  large 
damping  ratios.  Equation  (l6)  can  be  written  then  as 

(TqS  f  l)(TqS  +  l)(Tb2sa  +  2^8  +  1) 

4°  “  ^  S2(T*S  +  1)(TXS  +  1)[(TXS  +  l)  (Ty2s2  +  2£yTyS  +  l)  ] 

(18) 

It  is  now  possible  to  choose  the  parameters  in  equation  (17)  to  match  those 
of  equation  (l8);  the  specific  values  are  given  in  table  III. 

The  remainder  of  the  system  consists  of  the  tracking  radar  whose 
response  should  be  made  relatively  fast  for  this  design  approach.  Further, 
from  the  plots  presented  in  figures  8  and  9  it  is  apparent  that  the  track¬ 
ing  radar  should  be  designed  so  as  not  to  alter  appreciably  the  spectral 
distribution  of  error  or  control  motions.  Thus  a  break  point  or  natural 
frequency  of  several  cycles  per  second  would  be  satisfactory. 

Comparisons  and  other  considerations.-  In  a  REAC  simulation  of  this 
system  its  operation  was  examined  in  further  detail.  Of  particular  impor¬ 
tance  in  this  examination  was  the  control  deflection  and  the  associated 
servo  energy.  It  was  found  that  the  control  deflections  were  effectively 
restrained  to  within  the  linear  range  and  as  a  result  the  servo  energy  was 
greatly  reduced.  For  a  given  servo  the  power  expended  is  proportional  to 
the  time  average  of  the  sum  of  the  absolute  displacements  of  control 
motion.  Thus  the  average  servo  power  over  a  sufficient  time  interval  of 
duration  t  can  be  determined  from 


average  servo  power  -  — ^ — 1 

For  this  system  there  is  more  than  a  two-thirds  reduction  in  required  servo 
energy  over  that  for  the  system  with  smallest  servo  energy  requirement  dis¬ 
cussed  in  reference  3.  The  saving  in  servo  energy  becomes  even  greater 
when  compared  to  systems  for  which  the  error  performance  becomes  progres¬ 
sively  worse  than  the  optimum.  The  control-deflection  time  history 
obtained  from  the  REAC  was  used  to  show  that  the  required  control  rates 
are  easily  attainable  since  they  rarely  exceeded  3  to  3*5  radians/sec. 

This  system  has  other  virtues.  For  example,  there  is  also  a  reduction 
of  voltages  within  the  circuit  to  reasonable  and  easily  obtainable  values. 
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From  this  simulation  work,  then,  two  facts  are  apparent.  First,  It 
Is  clear  that  the  validity  of  applying  only  a  control-motion  restriction 
In  the  analysis  hac  been  verified,  since  the  control  rate  and  system  vol¬ 
tages  did  not  saturate.  Second,  It  is  apparent  that  in  comparison  with 
systems  of  reference  3  in  which  limiting  occurred,  the  present  system  not 
only  produced  a  smaller  error  hut  was  accompanied  by  other  desirable 
effects  such  as  decreased  servo  energy  and  circuit  voltage  requirements. 

The  missile-control  system  discussed  In  the  above  paragraphs  is 
Intended  only  as  an  exsmple  designed  to  achieve  optimum  performance 
against  noise.  It  is  clear  that  since  the  system  operation  is  essen¬ 
tially  linear,  all  designs  which  approximate  the  desired  Yo  would  give 
approximately  the  same  performance.  However,  it  is  often  necessary  to 
satisfy  other  requirements  which  are  related  to  the  specific  design  objec¬ 
tives  of  the  system  and  which  may  favor  certain  designs  or  system  con¬ 
figurations.  Because  of  the  freedom  in  choice  in  distributing  the  opti¬ 
mum  filtering  in  the  system,  the  optimum  characteristics  specified  by  the 
theory  do  not,  in  general,  place  an  inherent  limitation  on  satisfying  other 
requirements.  For  example,  requirements  on  transient  response  during  the 
beam  entry  phase  may  dictate  a  design  in  which  more  of  the  optimum  filter¬ 
ing  is  located  in  the  tracking  radar.  Similarly,  requirements  of  simpli¬ 
city  on  certain  parts  of  the  system  may  dictate  yet  a  different  design. 
Because  these  and  similar  requirements  are  related  to  the  detailed  design 
objectives,  an  investigation  of  such  factors  is  beyond  the  scope  of  this 
report. 


CONCLUDING  REMARKS 


This  study  has  considered  the  application  of  Newton's  modification 
of  the  Wiener  filter  theory  to  the  choice  of  optimum  transfer  function 
for  the  beam-rider  guidance  system.  The  analysis  has  been  applied  to  a 
typical  variable-incidence  missile  at  a  given  flight  condition  attacking 
a  maneuvering  target  in  the  presence  of  glint  noise.  By  minimizing  the 
mean-square  error  with  a  side  restriction  on  the  mean -square  value  of 
certain  quantities,  limiting  effects,  of  primary  importance  in  missile 
guidance,  are  largely  eliminated. 

Although  many  limiting-type  nonlinearities  are  of  importance,  this 
study  has  shown  that  the  critical  nonlinearity  is  due  to  control-deflection 
limiting,  since  a  restriction  applied  to  the  deflection  sufficient  to  con¬ 
fine  control  motions  to  within  the  linear  range  also  satisfactorily  pre¬ 
vented  other  types  of  limiting.  Results  have  been  given  which  illustrate 
the  effect  of  the  control-motion  restriction  on  the  minimum  rms  error. 

Zero  restriction  corresponds  to  the  Wiener  theory  whose  results  are 
unattainable  since  impossibly  large  control  motions  are  demanded.  How¬ 
ever,  as  the  surface  deflection  is  reduced  by  means  of  an  increasing 
restriction,  the  minimum  error  increases  slowly  until  at  a  realistic 
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value  of  the  control  deflection  the  error  has  increased  only  a  few  feet. 
Thus  it  is  clear  that  the  optimum  performance  predicted  by  the  Wiener 
theory  can  still  be  approached  closely  by  confining  operation  of  the  sys¬ 
tem  to  a  realistic  linear  range. 

The  transfer  functions  of  the  optimum  guidance  system  have  been  given 
and  have  been  compared  to  the  corresponding  transfer  functions  from  the 
Wiener  theory.  Differences  vere  shown  to  be  the  result  of  added  filtering 
terms  demanded  by  the  restriction  on  control  motion.  The  various  filtering 
factors  can  be  traced  to  their  origin.  It  was  shown  that  terms  due  to 
basic  Wiener  theory  appear  in  the  solution,  and  that  the  additional  terms 
can  be  associated  with  the  control-motion  restriction  and  the  missile  aero¬ 
dynamics. 

For  the  optimum  system,  performance  against  variable  noise  magnitude 
was  given  and  was  compared  to  Wiener  theory.  As  an  example  of  one  possible 
design  which  achieves  these  results,  a  missile-control  system  was  synthe¬ 
sized  to  match  the  optimum  transfer  function.  In  a  REAC  simulation  the 
operation  of  this  sytem  was  examined  in  further  detail.  It  was  found  that 
the  magnitude  of  circuit  voltages,  and  control-motion  rates  are  greatly 
reduced  in  comparison  with  systems  in  which  limiting  occurs.  Similarly, 
the  required  servo  energies  are  appreciably  reduced,  resulting  in  about 
a  two-thirds  saving. 
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APPENDIX  A 


DETERMINATION  OF  THE  OPTIMUM  MODIFIED  WIENER  FILTER 


The  solution  for  the  optimum  compensating  network  Hco  illustrated 
in  figure  3  is  given  by  equation  (6)  in  the  text  as  follows: 


where 


Hco(iw)  *  r  •**'  r  -  « 

ncoU  '  2*A+(u>)  J0  A' (a) 

A(w)  =  ^Hf  (iai)Hf  (iu>)  +  pj«ii(t»)) 


>  (Al) 


In  this  equation  4>ii  is  defined  as 


®ii  ■  ®T  +  ®N  (A2) 

where  4^  and  4>N  represent  the  spectral  density  of  the  target  motion  and 
noise |  respectively.  The  quantities  A+  and  A~  are  defined  as  the  factors 
of  A  with  poles  and  zeros  in  the  upper  and  lower  half -planes,  respec¬ 
tively.  Thus 


A  =  A+A" 


(A3) 


For  the  case  to  be  considered  here  the  target  motion  and  noise  will  be 
taken  as 


kaiji2 

+  k2) 


(a4) 


4>n  -  N 


(A5) 


It  will  be  noted  that  the  noise  spectrum  has  been  approximated  here  by  a 
constant  in  order  to  reduce  the  complexity  of  the  calculations.  As  shown 
in  reference  3,  this  approximation  is  valid  because  the  noise  defined  by 
equation  (l)  is  essentially  flat  compared  to  the  passband  of  the  optimum 
system.  It  should  also  be  pointed  out  that  use  of  equation  (A4)  leads  to 
certain  mathematical  difficulties  which  occur  in  the  process  of  evaluating 
the  right-hand  side  of  equation  (Al),  because  the  theory  requires  that  the 
poles  of  4^  not  be  located  on  the  real  axis.  To  avoid  these  difficulties 
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it  is  necessary  to  modify  the  target ^motion  spectral  density  to  the 
following 

kaT2 

*T  "  +  n1e)(<J2  +  ti22)(^  +  k2)  (A6) 

where  tj1  and  t\z  are  any  small  real  numbers .  The  solution  to  equation 
(Al)  will  then  be  a  function  of  lu,  tjx ,  and  The  desired  answer  is 
obtained  by  taking  lim  Y0  ,tj2)  as  and  t|2  -»  0.  However,  it  can 

be  shown  (see  ref.  3)  that  the  same  answer  can  be  obtained  more  simply  by 
taking  tj1  and  n2  equal  to  zero  as  would  be  obtained  from  the  rigorous 
process  described  above. 

Now  by  combining  equations  (A2) ,  (kk) ,  and  (A5) , 


®ii(a)  =  «T(a)  +  ®N(a) 

ag  +  k2a4  +  (ks«p2/itN) 
a4 (a2  +  k2) 


(a  -  api) (a  -  apg)(a  -  aps) (a  +  dpi) (a  +  ap2) (a  +  dps) 
a4(a  +  ik)(a  -  ik) 


N  [p(a)3[-p(-a)] 
a4 (a  +  ik)(a  -  ik) 


(A7) 


where  the  apm  represent  roots  in  the  upper  half -plane.  The  polynomial 
p(a)  is 

3 

p(a)  =  II  (a  -  cO 
m**i  ^ 


s  a3  +  b2a2  +  bxa  +  bo 


(A8) 


where 


b2  =  -(api  +  Op2  +  tips) 

bi  =  a.p1ap2  +  apxaps  +  ap2aps  • 
bo  =  -cipxcip2apg 


(A9) 
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In  a  similar  manner 


-p(-a)  =  a3  -  b2a2  +  bia  -  bo  (A10) 

An  expression  to  be  used  later  can  be  obtained  from  these  preceding 
equations  as  follows 


a8  +  k2a4  +  "“jjj-*  (a.3  +  b2a2  +  bia  +  bo)  (a3  -  b2a2  +  bxa  -  bo) 
from  which 

-bo2  =  ISt 


nN 


(All) 


The  quantity  Hf  in  equation  (Al)  represents  the  fixed  network  or 
aerodynamics  and  its  general  form  (ref.  7)  is  given  by 

Hf(la)  .  3«  .  1  -  T^a2)  (A12) 

HfU  J  6  Ti5  a2[(l  -  Ta2a2)  +  i(25aTaa)]  K  ’ 


Hr(ia)  =  -  -i*  •  I, \  -  -  ^Tkai] .  (AI3) 

UfUa;  a2[  (l  -  Ta2az)  -  i(2^aTaa)] 


and 


Thus 


=-rr^r,  U  \  1  [Tb4a4  +  (^^gTbg  -  gTbg)ag  4  1] 

Hf(ia)Hf(ia)  -  ^  ^^4  +  (^2^2  .  2Ta2)a2  +  1] 

Accordingly  it  can  be  shown  that  in  equation  (Al) 


(A14) 


Hf  (ia)Hf  (ia)  +  p 


pTs4Ta4aa+pTs4(4gagTag-gTag)aa^(Tb4+pTs4)a4+(4^2Tb2-2Tb2)ag+l 
Ts4w4[Ta4a4+( 4£a2Ta2 -2Ta2 ) a2+l ] 

4  4 

II  (a-aqn)  II  (a+aqn) 


=  dT  4  _ 2=1 _ BSi - 

P  a  a4 [ Ta4a4+ ( 4  ?a2Ta2 -2Ta2 ) az+l ] 


&  pTa4  —t 


q(a)q(-q) 


a4 [ Ta4a4+ ( 4  ?a2Ta2 -2Ta2 ) a2+l ] 


(A15) 
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This  expression  can  be  split  Into  two  factors  with  poles  and  zeros  In  the 
upper  and  lower  half -planes .  From  the  preceding  definitions  q(a)  and 
p(a)  both  have  their  zeros  In  the  upper  half -plane.  Thus 


.+,  v  _4_ _ q(a)p(a) _ 

n  a  =  P  a  J  a4(l  -  Ta2a2  +  12?aTaa) (a  -  lk) 


(A21) 


A"  (a) 


q(-a)[-p(-a)] 


a4(l  -  Ta2a2  -  12taTaa) (a  +  lk) 


(A22) 


Now.  let  us  consider  the  evaluation  of  the  first  integral  in  equa¬ 
tion  (Al) .  This  Integral  can  be  evaluated  by  combining  equations  (A4) , 
(A13) ,  and  (A22)  as  follows: 

IX  s  r  ^.(a)e^  da 

j-  00  A' (a) 

kaT2  r”  (1  -  Tb2a2  -  12^a)elat 
Tb2*  J.« o  a2(ot  -  ik)q(-a)  [  -p(-a)  ] 


T 7* 


/ 


(1  -  T^a2  - 


ya  " 12^ 


a)elat 


da 


r 


Tg  It  00 


a2(a  -  ik)  II  (a  +  aqn)  II  (a  +  a™,) 
n=i  m=i 


f (a) da 


(A23) 


where  f(a)  is  defined  as  the  Integrand.  By  considering  a  to  be  a  com¬ 
plex  variable,  the  integral  can  be  evaluated  by  a  contour  integration  indi¬ 
cated  by  the  following  sketch: 
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Thus  it  can  he  shown  that 

21kaT2  (7i  +  72)boc0(iu))2  +  [krtbpCp  +  (l/k)]lfa>  +  1 
2  "  Ts^boCo  (iw)2(k  +  iw) 

21kaTg  T^do))2  +  2^T„(i<»)  +  1 

S  Ts^boco  (ito)2(k  +  iu>) 


It  should  be  noted  that  since  the  q's  were  assumed  to  be  zero,  questions 
concerning  the  existence  of  the  above  integral  arise.  However,  as  indi¬ 
cated  in  reference  3,  when  the  analysis  is  made  without  this  assumption, 
the  integral  in  equation  (A27)  does  exist,  and  the  limit  of  this  integral 
as  the  q*s  approach  zero  becomes  precisely  equation  (A28) . 

The  coefficients  can  be  expressed  in  a  more  convenient  form  by  elimi¬ 
nating  the  intermediate  parameters  y1  and  yz  by  means  of  definitions 
given  in  equations  (A24)  and  (A25) .  The  following  results  are  then 
obtained: 

2£aTa  =  k7jb0c0  + 


-  -2£bTb 


^  bpc^  +  bxCp 

b0co 


(A29) 


Also 

TaZ  -  (7x  +  72)b0c0 

b0c0(Thgk2  +  2ShThK  *  1)  _  kbpc^^Tb  +  ik(boCl  +  bxc0)  _  1_ 
k2q( -ik) [ -p( -ik) ]  k^boco  k2 

(A30) 

For  the  missile  parameters  of  table  I  and  the  value  of  k  given  on  page  28 

Tb2k2  +  2W  +  1-1 


so  that  very  nearly 


Ta2 


k2q(-ik)[-p(-ik)] 


1_ 

k2 


(A31) 
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Nov,  from  equation  (Al)  the  optimum  compensating  network  can  he 
determined  as  follows: 

H=°(1“>  ■  dfo  <*3a) 

Substituting  equations  (A2l)  and  (A28)  Into  (A32)  gives 
Hco(iw) 


kar2  ( la)) 2  [  T„2  ( lui)  2-tfi  Cn.Tcr  ( 1(»)  +1  ]  [  Tr2  ( lfa>) 2 +2  kTa  ( Ifai)  +1] 

aNpTs2Ta4hoco  p(u))q(w) 


(A33) 


To  reduce  this  expression  to  a  more  standard  form  It  can  be  seen  from 
equation  (Al6)  that 

4 

q(u)  -  (w  -  aqn) 

"  Co  nt  ("5^  +  X)  (A34) 

Similarly  from  equation  (A 8) 

3 

p(w)  =  nSi  (w  -  ttpm) 

(A35) 

Also  by  using  the  definitions  given  in  equations  (All)  and  (A19)  the  com¬ 
pensating  network  can  be  shown  to  reduce  to 


HcoUu)  -  Ts2  <lu>8t W“)a  +  yaTq(ito)  +  HtVd^)2  *  2^*1**)  +  H 

(A36) 

By  the  substitution  of  the  conventional  complex  frequency  s  for  lea 


Hco(s)  =  Ts* 


;(Tn2s2  +  2t '  Twb  4  l)(Ta2s2  +  2^,Ths  +  l) 


(A37) 
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Because  and  (;x  are  close  enough  to  unity  (table  II)  it  is  possible 
to  factor  the  corresponding  terms  approximately  to  give 

„  ^♦l)(TBm^l)(Vs«*aebV..H) 

^  S2(T^B  +  l)(TxS  +  l)  [  (TXS  +  l)(Ty2Sa  +  2£yTyS  +  l)  ] 

(A42) 

The  numerical  evaluation  of  the  optimum  system  has  been  carried  out 
for  the  following  values  of  target  motion  and  noise: 

k  =  0.4  flip  «  1  g 

T  =  5  sec  N  =  15  ft2/radian/sec 

From  this  example  p  will  be  taken  to  be  104  since,  as  shown  in  the  text, 
this  value  results  in  a  near-optimum  restriction  on  control  deflection. 

Evaluation  of  the  numerator  in  equation  (A7)  gives 

ka_,2 

a®  +  k2a4  +  ^  =  a®  +  0.l6  a4  +  8.8 
3  3 

=  (a  -  <W  m31  (a  +  apn) 

where 

apl  =  i  1.456 

apa  =  1.427  erp(i  0.543) 

aps  m  -1.427  exp(-i  O.543) 

which  are  all  located  in  the  upper  half  of  the  complex  a  plane.  Then 
from  equation  (A9) 


b2  =  -i  2.94 
hi  «  -4.201 
b0  =  i  2.97 
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Similarly  from  equation  (A15)  for  p  ■  104, 
pT84Ta4a8+pT84(4Ca2Ta2-2Ta2)a®+(Tb4+pTs4)a4+(45b2Tb2-2Tba)a2+l 

=  (22.58xl0-e)a8-(7.475xl0-H)a®+(6.267X10-*)a4-(0.006094)a2+l 

4  4 

*  pTs^a4  ^  (a-ctqn)  ^  (o+oqn) 


Solving  for  the  roots  gives 

otqi  =  3-593  exp(i  0.8o4) 
a<32  =  -3-593  exp(-i  0.804) 
otqg  =  12.92  exp(i  0.058) 
aq4  =  -12.92  exp(-i  0.058) 

which  again  are  all  in  the  upper-half  a  plane.  From  the  definitions  in 
equation  (A17) , 


cs  -  -i  6.677 
c2  =  -187.61 
cx  =  i  883.40 
co  -  2155 

Now  from  the  values  of  the  roots  apm  and  aqn  as  well  as  the  h  and  c 
coefficients,  the  optimum  compensating  network  can  he  calculated  by  means 
of  equation  (A38) .  For  this  purpose  it  is  necessary  to  evaluate  the  left- 
hand  sides  of  equations  (A29) ,  (A31) ,  and  the  parameters  following  equa¬ 
tion  (A38) .  Results  of  these  calculations  are  given  in  table  II  in  terms 
of  damping  ratios  and  undamped  natural  periods.  Substitution  of  these 
parameters  in  equations  (A38)  and  (A40)  gives 
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APPENDIX  B 

MISSILE-CONTROL -SYSTB*  APPROXIMATION  TO 
THE  OPTIMUM  TRANSFER  FUNCTION 


The  optimum  linear  system  is  illustrated  in  figure  11.  In  the  follow¬ 
ing  section  the  system  equations  are  derived. 

Using  figure  11  it  can  he  shown  that 

J_  m  _ kfl(^s  +  l)(Tss  +  1) _ 

%  T^CTjS  +  l)(Ta2s2  +  2?aTas  +  l)(Tas  +  l)  +  kBkss(TMs  +  l)(T2s  +  l) 

(Bl) 

To  simplify  equation  (Bl)  it  is  convenient  to  make  T3  *  %.  This  assump¬ 
tion  is  not  essential  hut  its  use  leads  to  simpler  equations.  In  certain 
cases  where  gust  disturbances  are  serious  it  may  he  more  desirable  to 
choose  T3  small.  With  the  former  choice 

-1 _ _ _ _ 

%  ■[(TaTaaT1)8*+(TdTa«4®d2{aTaT1)8#+(TdT14Td25aTa*Bk3Ta)s+(Td4kBk3)] 


k8 

Td+ksk3 


(T^S+l) 

s(ass+bs2+cs+l) 


(B2) 


where 


Td  +  kBka 

=  Td(Ta2  f  2gaTaTt) 

Td  +  ksk3 

„  ^,k.T-a.%.  jL'EiTdl,  1  hfca?.z 

Td  +  ksk3 


(B3) 


,3-1 y 
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Thus  the  entire  open-loop  transfer  function  can  be  written  as 

534  k2ksV  (T10s  +  1)(T4b  +  l)(Tb2s2  +  2^8  +  l) 


{  Td  +  kska  b® (TjjB  +  1)(T5s  +  l) (as®  +  bs2  +  cs  +  1) 


By  comparison  with  equation  (A42)  it  can  be  seen  to  be  of  the  same  form  as 
the  optimum  transfer  function.  Thus  the  following  correspondences  between 
quantities  in  equations  (B4)  and  (A42)  are  appropriate: 


Tn  «->  T* 


Also  it  is  apparent  that 


Tio  «-*  Ta 


T4  *->  Ta 


TxTy2 


-Tyz  +  Tx2?yTy 
’  Tx  +  2?yTy 


For  these  values  of  a,  b,  and  c,  equations  (B3)  can  be  solved  for  the 
remaining  parameters  to  give  . 

Tl  =  — 

11  bTa  -  a2?a 


ksk3  = 


TdTa2Ti  -  aTd 


cTa  +  ckBka  -  2;aTftT^  -  T,  T, 


/ 
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The  value  of  k2ks  Is  determined  from  the  desired  gain: 

kgks  <->  ^  (b8) 

Since  only  k2ks  and  kska  are  specified,  one  of  the  three  gains  k2,  ks, 
or  ks  may  he  chosen  arbitrarily,  subject  to  the  condition  that  voltage 
limiting  does  not  occur.  From  the  above  equations  the  system  parameters 
have  been  determined  and  are  tabulated  in  table  III. 
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TABLE  I.-  SUMMARY  OF  AERODYNAMIC  PARAMETERS  FOR  EXAMPLE  MISSILE 


Parameter 

Value 

Ta 

0.0775 

Tb 

.0552 

Td 

2.087 

Ts2 

.00079H 

.846 

£a 

.0536 

.0220 

TABLE  II.-  SUMMARY  OF  PARAMETERS  OF  OPTIMUM  TRANSFER  FUNCTIONS 


Parameter 


Modified  Wiener 
theory  value 

Wiener  theory 
value 

1.13 

0.925 

.805 

.765 

.687 

.687 

.700 

.700 

•  519 

.519 

2.5 

2-5 

3-2 

7.42 

.278 

— 

.720 

— 

.0774 

— 

.0581 

— 

.184 

— 

.722 

.0787 

— 

.0513 

— 
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TABLE  III.-  SUMMARY  OF  PARAMETERS  OF  OPTIMUM  SYSTEM 


Parameter 

Value 

k2 

0.0538 

T4 

1.13 

T10 

1.13 

Til 

2-5 

t5 

.184 

k8 

.0463 

Ti 

.186 

*3 

-.674 

t2 

•  373 

Ts 

i  .846 

Figure  2.-  Block  diagram  of  general  filter  problem  with  restrictions 
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Figure  4.-  The  beam-rider  system. 


Figure  7*-  Concluded. 


Frequency  ,cps  . 

Figure  9-“  Spectral  density  of  error  components 


